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Compact U(l) gauge theories in 2 + 1 dimensions and the physics 
of low dimensional insulating materials 
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Abstract. Compact abelian gauge theories in d = 2 + 1 dimensions arise often as an effective field-theoretic 
description of models of quantum insulators. In this paper we review some recent results about the compact 
abelian Higgs model in d — 2 + 1 in that context. 

PACS. 11. 15. Ha Lattice gauge theory - ll.10.Kk Field theories in dimensions other than four 



1 Introduction 

Effective abelian gauge theories of doped Mott insulators 
are very popular in the condensed matter literature. De- 
spite the enormous amount of papers published over the 
last ten years, there are still many points of disagreement 
among researchers, mostly related to poorly understood is- 
sues [T1I21I31I3] . One key point is the interplay between the 
confinement/deconfinement properties of the many the- 
ories available, and spin-charge separation. It is usually 
believed that if deconfinement occurs, then spin-charge 
separation would also occur. However, it remains a con- 
troversial issue if deconfinement is possible for gauge fields 
coupled to matter with the fundamental charge. The prob- 
lem here is that most U(l) effective gauge theories are 
studied in d = 2 + 1 space-time dimensions. For this di- 
mensionality t/(l) gauge fields are strongly confining. For 
instance, if we neglect the coupling to matter fields we ob- 
tain that pure compact Maxwell theory permanently con- 
fines electric charges in d = 2 + 1 dimensions 5j . This is in 
contrast with pure compact Maxwell theory in d = 3 + 1 
dimensions where a deconfinement transition occurs [jy. 
A natural question to ask in this context is whether the 
coupling to matter fields is able to change the permanent 
confinement in 2 + 1 dimensions. One example where the 
answer to this question is affirmative corresponds to the 
case where bosons are coupled to an abelian gauge field 
with integer charge q > 1 |Zj . This is most easily under- 
stood by considering the lattice gauge theory version of 
such a model, whose action is given by 

S = -/35^[cos(V„fl 4 - qA itl ) -1]-kJ2 [cos(iV) - 1], 



where we have V^flj = 



H+fj, 



a - 0i and F iliU = V M A it/ 



(1) 



X7 v Ain. The above action corresponds to the compact lat- 



tice abelian Higgs (CLAH) model. When j3 — ► oo the the- 
ory becomes a Z q gauge theory which is known to have a 
deconfining transition in 2 + 1 dimensions. For k — > oo we 
obtain the XY model, which also undergoes a phase tran- 
sition. The critical points corresponding to these two lim- 
iting cases are connected by a critical line separating the 
confining from the deconfining phase. There is no Coulomb 
phase in 2 + 1 dimensions. The phase diagram and the 
critical properties of the action were studied in de- 
tail using large scale Monte Carlo simulations in Ref. 
The CLAH model appears in many contexts in condensed 
matter physics. We shall cite only two important recent 
examples. In the q = 2 version it arises as an effective the- 
ory for the two-dimensional quantum Heisenberg antifer- 
romagnet (QHA) !) . There the gauge field originates from 
the Berry phase. Another recent example corresponds to 
the strongly correlated limit of a bosonic model considered 
by Senthil and Motrunich JOj- In this context the CLAH 
model describes a transition from an ordinary Mott insu- 
lating phase to a fractionalized insulating phase 

However, bosonic actions like the one in Eq. are 
only part of the complete effective action associated with 
a doped Mott insulator. It usually contains also fermions 
coupled to a gauge field. [TTUT^j The traditional route 
starts with a slave boson representation of the t — J model 
where the projected electron operator is written as a com- 
posite particle, Ci a = &j/j a , where bi is an auxiliary bo- 
son and fi a is an auxiliary fermion. The auxiliary fields 
obey the constraint b\bi + ^2 a fj a fi a = 1 at each lattice 
site. In the resulting effective action the auxiliary fields 
interact through a compact gauge field. The most popular 
way of doing things attribute the charge of electron to the 
bosonic field. Thus, the fermion would have no charge and 
would only carry the spin degree of freedom. If the effec- 
tive model undergoes a deconfinement transition fermions 
and bosons will have an independent dynamics and in this 
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way we can say the spin and charge are separated. Nayak 
pointed out recently that things are not so simple be- 
cause the assignment of the whole electron charge to the 
auxiliary boson is completely arbitrary. We could well say 
that the boson carries charge Se while the fermion carries 
charge (1 — 6)e. The constraints of the theory, which are 
enforced by the coupling to the gauge field, would then en- 
sure that the physics is unaffected by this arbitrary choice. 
This scenario implies that the auxiliary fields should not 
be associated with the physical charge and spin excita- 
tions of the theory. The auxiliary fields introduced this 
way are not part of the physical spectrum. This situation 
is reminiscent of the analysis made by Mudry and Fradkin 
some time ago. [13] 

In this paper, we review some recents results on the 
model (T1J. We start discussing in the Section [21 the finite- 
size scaling (FSS) analysis of the model in the q > 1 case. 
We employ a new and very accurate method for extracting 
the critical exponents. In Scction[3]we discuss the phase di- 
agram, which is obtained for q = 2,3, 4, 5. In Section0]we 
discuss the q = 1 case, whose phase diagram is presently 
a matter of intensive debate [r4lll5lll6| . 

2 The universality class of the deconfinement 
transition for q > 1 

In 1981 Bhanot and Freedman FT7] made a finite-size scal- 
ing (FSS) analysis of the CLAH model for q — 2 and 
d = 2 + 1. They obtained the phase diagram of the model 
but the critical exponents were not calculated. In order 
to determine the universality class of the model it is nec- 
essary to compute the critical exponents. The accurate 
determination of the critical behavior is not an easy task. 
In principle, on the basis of the results of Fradkin and 
Shenker [7j, we might think that the universality class 
will be the same as of the Z q spin model, except for the 
limit k — ► oo where the universality class is obviously that 
of the XY model. If we assume that this is indeed so, we 
obtain that the critical exponents for the q = 2 case have 
Ising values. However, we have recently shown through 
large scale Monte Carlo simulations that the situation is 
more complicated [S]. For example, for q — 2 we have ob- 
tained that the exponents are those of the Ising model in 
a large part of the phase diagram in the K/3-plane, but 
there is a finite interval k\ < k < K2 where the exponents 
vary continuously, before reaching XY values for k > k%. 
Therefore, it seems that the model features a fixed line 
rather than a fixed point and belongs to a new universal- 
ity class. 

The FSS analysis of Bhanot and Freedman fT2] relies 
on the second moment of the free energy, i.e., the specific 
heat. In this case it is well known that very large system 
sizes are needed in order to identify a genuine critical be- 
havior to high accuracy. We have shown that much better 
results are obtained by performing a FSS analysis based 
on the third moment of the free energy. Such a procedure 
is not only more accurate: it allows also for an independent 
determination of the exponents v and a, providing in this 




Coupling constant 



Fig. 1. Generic third moment M3, showing how FSS is used 
to extract a and v 

way a check of hyperscaling |Sj . The reason for this lies in 
the double-peak structure of the third moment. The FSS 
in this case is such that the width between the negative 
and positive peak scales as L~ x l v , while the height of the 
positive peak with respect to the negative one scales as 
£(i+a)/r/_ xhis scaling behavior is shown schematically for 
a generic third moment M3 in Fig. F2J 

The critical exponents as a function of k/@ are shown 
in Fig.EJfor q = 2. It is clearly seen that the Ising critical 
behavior connects to the XY behavior through a regime 
of continuously varying critical exponents. 



3 The phase diagram for q > 1 

We have obtained the phase diagram for q — 2,3,4,5. 
Among these values of q, only the q = 3 case exhibits 
a first-order phase transition up to some point where it 
changes to second-order. This behavior is consistent with 
the fact that the three-state Potts model exhibits a first- 
order transition in three dimensions. This model is dual 
to the Z 3 gauge theory. The phase diagram for q = 3 case 
therefore features a tricritical point. The tricritical point 
in the phase diagram is estimated as (/3tri, ttri) = (1.23 ± 
0.03, 1.73 ± 0.03), corresponding to a ratio Ktrj/Ari = 
1.39 ± 0.06. At this point it is worth making the follow- 
ing remark. The non-compact version of the model iTjJ is 
known to exhibit a second-order phase transition for all 
values of q. As matter of fact, the value of q is not im- 
portant in that case, and can be absorbed into a redefini- 
tion of the gauge field. The universality class corresponds 
to the so called "inverted" XY transition |TH|. However, 
if the radial part of the scalar field is allowed to fluctu- 
ate, thus generalizing the non-compact version of (TJ3), then 
it is possible to show using duality arguments that the 
resulting model has a tricritical point approximately at 
Ktri//3tri ~ 0-64 This estimate was confirmed recently 
by large scale Monte Carlo simulations [20| . In the case of 
the q = 3 CLAH model there is no need to consider the 
fluctuations of the radial component of the field to obtain a 
tricritical point. The origin of this tricritical point is com- 
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Fig. 2. a) (1 + a)/v from FSS finite-size of M 3 for Eq. JTJ for 
g = 2. Note the variation relative to the ^2-limit 1.76 (dotted 
horizontal line) and the (7(l)-limit 1.467 (solid horizontal line). 

b) Same for the exponent v, computed directly from M3 (A) 
and combining results for (1 + a)jv with hyperscaling (O)- 

c) a as computed directly from M3 (A) and using results for 
(l + a)/^ with hyperscaling (O)- The maximum and minimum 
in a) have been obtained by crossing the critical line along the 
trajectory /3(k) = f3 c + a (ft — k c ) with a — 00 (A), a = 1 
(■), and a — —1 (a) using /3 C — 0.665, k c = 2.125 (max.), and 
0c = 0.525, k c = 5.0 (min.). 



pletely different from the one of the non-compact model. 
Anyway, it is interesting to note that the ratio Ktri//3tri in 
the q = 3 CLAH model is more than twice larger than the 
corresponding ratio in the non-compact case. 

In the context of the model discussed recently by Senthil 
and Motrunich the phase diagram for the q = 2 

case has the following physical meaning. The confined 
phase corresponds to gap to integer charge excitations 



and is therefore interpreted as an ordinary Mott insulat- 
ing phase. The deconfined phase, on the other hand, is the 
same as the Higgs phase. This gives a gap to half-integer 
charge excitations and is interpreted as a fractionalized 
insulator. The model can in principle be artificially real- 
ized with present day technology by building an array of 
Josephson junctions of a particular type jl()| . 




Fig. 3. The phase diagram for the d = 3 compact abelian 
Higgs model in three dimensions for q = 2, 3, 4, 5. All lines are 
critical for all values of k except for the case q = 3, which is 
first order for k < ft tr i and second order otherwise. The thick 
solid portion of the q = 3 line is a first-order transition line. 
The critical line approaches the 3DXY value /3 C = 0.453 as 
k — > 00 for all integer values of q > 1. 



4 The q = 1 case 

The q — 1 case was generally believed to not exhibit any 
phase transition at all. The coupling to matter fields would 
not change the scenario of the d = 2 + 1 compact Maxwell 
theory, where it is known that no phase transition takes 
place the electric charges would never deconfine. A 
research initiated with the study of the QHA showed that 
things seem not to be so simple 1 1 . The main point is that 
by integrating out the matter fields we obtain to quadratic 
order in the gauge field fluctuations the effective action 



S e g (X / d xF, 



1 



--F„ 



(2) 



Therefore, while in the absence of matter fields the monopoles 
interact through a 1/r potential, in the presence of matter 
fields this behavior changes to — lnr. In the former case 
it is well known that upon dualizing the theory a sine- 
Gordon action is obtained. In the latter case, however, an 
anomalous sine-Gordon action arises |14| : 



S = 



1 



8tt 2 K 



d 3 x[ip{-d 2 ) 3/2 ip-2zcostp], (3) 
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where K = 1/g 2 , with g being the gauge coupling. Thus, 
we arrive at a behavior similar to the two-dimensional 
case. By simply looking at the scaling behavior of the fu- 
gacity of the monopole gas, one would be lead to conclude 
that a Kosterlitz-Thouless (KT) phase transition happens 
in three space-time dimensions. Note that this KT tran- 
sition is at zero temperature and has nothing to do with 
other KT transitions obtained by dimensional reduction 
due to temperature effects. However, it is not enough to 
study the scaling of just the fugacity. Since we are in 
d — 2 + 1 a more thorough analysis is needed in order 
to fully establish the KT behavior. To this end, it is nec- 
essary to also know the scaling behavior of the "stiffness" 
K, since the scaling of the fugacity and the scaling of the 
stiffness mutually influence each other in a subtle way. 

For the usual sine-Gordon model in d dimensions we 
obtain the following coupled recursion relations for the 
fugacity and the stiffness [2] 



dK- 



dl 
dy 



4ttV - (2 - d)K~ 



^=[d-2n 2 f(d)K] y, 



(4) 



(5) 



where f{d) = (d - 2)r[(d - 2)/2]/(47r) d / 2 . For d = 2 
the above equations imply the existence of a line of fixed 
points characteristic of the KT transition. For d = 3, how- 
ever, no fixed point is found and therefore no phase transi- 
tion takes place. This is consistent with Polyakov's result 
for the compact Maxwell theory in d = 2 + 1 5 j . The issue 
here is how the above recursion relations are modified in 
the case of the anomalous sine-Gordon model, Eq. J3J ■ To 
investigate this in great generality we will consider screen- 
ing in a Coulomb gas of monopoles with a propagator of 
the form in 2 < d < 4 dimensions with a > 0. 

The particular case d = a = 3 corresponds to our anoma- 
lous sine-Gordon theory. The bare potential is given by 
U (r) = -An 2 KV(r), where 



V(r) = 



r 



' d-c 
v 2 



-[{At) 



1], 



(6) 



2< J ir d / 2 r(a/2 

with A being an ultraviolet cutoff. Note that for d = a we 
obtain a potential oc In r. The bare electric field is given by 
E (r) = -A-K 2 KA{d,a)r' J - d - 1 /r%~ d , where r Q = 1/A and 
A(d,a) = {d - a)r[{d - a)/2]/[2 a Ti d / 2 r{a/2)]. The bare 
electric field is renormalized by the other dipolcs which are 
treated as a dielectric medium. The renormalized electric 
field is then given by 



E(r) 



4ir 2 K A(d,a) 



-d-i 



e(r) 



(7) 



where e(r) is the scale-dependent dielectric constant of 
the medium. This problem can be solved self-consistently 
for a renormalized potential U{r) derived from the above 
electric field |14| . In this way we obtain the recursion re- 
lations 



dK- 



dl 



4ttV 



(a - d)K~ 



(8) 



dy 
dl 



r, y - 2ir 2 A(d, a)K] y, 



(9) 



where the anomalous dimension of the fugacity is given 
by T] y = (ct — 2)/2. The above equations reduce to Eqs. (J2J 
and (0) when a = 2. The case relevant to the q = 1 CLAH 
model is a = d = 3. In this case Eqs. © and © are very 
similar to the usual KT recursion relations, except for the 
presence of the anomalous scaling dimension of the fugac- 
ity, rj y which is nonzero in our case and given by rj y = 1/2. 
Thus, we obtain that the monopole gas undergoes a KT- 
like phase transition. This is in contrast with Polyakov's 
compact Maxwell theory where the monopoles are always 
in the plasma phase. Using the usual duality arguments, 
we identify the plasma phase with the confined phase for 
the electric charges. The dielectric phase of the monopole 
gas is accordingly identified with the deconfined phase for 
the electric charges. 

At this point an important remark is in order. The 
above screening analysis was made in real space, in the 
spirit of the original Kosterlitz and Thouless paper |22j. 
Indeed, we have considered a local scale dependent dielec- 
tric constant. Screening arguments in momentum space 
usually involve a dielectric constant which is local in mo- 
mentum space, leading to an electric displacement vector 
of the form D(r) = J d d r'e(r — r')E(r'), while in our case 
we have simply D(r) = e(r)E(r). In the present context 
this may be understood by using a classical electrostatic 
argument. Let us consider the potential © in the case 
a = d, such that we have a In r-potential in d dimensions. 
Due to Gauss's theorem in d-dimensions, the field equa- 
tion is given by |22j V • (Er 2 ~ rf ) = Sdp(r), where Sd is the 
solid angle in d-dimensions. This can be cast in a more 
familiar form V • D = S^r^ -2 p(r) by introducing the "di- 
electric constant" £o(r) — (r/ro) 2 ~ d . Note that for d > 2 
this dielectric constant vanishes at large distances. This is 
rather unusual since in classical electromagnetism it can 
be shown that the dielectric constant is always greater 
than one. However, such a situation corresponds precisely 
to the anti-screening mechanism discussed many years ago 
in theories of confinement [21] • Thus, already at this level, 
even before taking the dipoles into account, we can define 
a kind of scale-dependent dielectric constant which is local 
in real space. 



5 Conclusions 

In this paper we discussed the deconfinement transition 
in the CLAH model, Eq. 0), in d = 2 + 1. To each value 
of the charge q corresponds a distinct universality class, 
in contrast to the non-compact model. While a decon- 
fined transition is expected for q > 1, the case q = 1 is 
still controversial [T^lll5lll6| and needs further study. One 
feature that also deserves further study is the possible ex- 
istence of a fixed line in the q = 2 case. The numerical 
results strongly suggest this possibility but multicritical 
fixed points are also possible. 
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